Noise-Induced Phase Space Transport in Two-Dimensional Hamiltonian Systems 
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First passage time experiments were used to explore the effects of low amplitude noise as a source 
of accelerated phase space diffusion in two-dimensional Hamiltonian systems, and these effects were 
then compared with the effects of periodic driving. The objective was to quantify and understand 
the manner in which "sticky" chaotic orbits that, in the absence of perturbations, are confined 
near regular islands for very long times, can become "unstuck" much more quickly when subjected 
to even very weak perturbations. For both noise and periodic driving, the typical escape time 
scales logarithmically with the amplitude of the perturbation. For white noise, the details seem 
unimportant: Additive and multiplicative noise typically have very similar effects, and the presence 
or absence of a friction related to the noise by a Fluctuation-Dissipation Theorem is also largely 
irrelevant. Allowing for colored noise can significantly decrease the efficacy of the perturbation, but 
only when the autocorrelation time, which vanishes for white noise, becomes so large that there is 
little power at frequencies comparable to the natural frequencies of the unperturbed orbit. Similarly, 
periodic driving is relatively inefficient when the driving frequency is not comparable to these natural 
frequencies. This suggests strongly that noise-induced extrinsic diffusion, like modulational diffusion 
associated with periodic driving, is a resonance phenomenon. The logarithmic dependence of the 
escape time on amplitude reflects the fact that the time required for perturbed and unperturbed 
orbits to diverge a given distance scales logarithmically in the amplitude of the perturbation. 



PACS number(s); 05.60.+W, 51.10.+y 05.40.+J 
I. MOTIVATION 

It is well known that a complex phase space containing 
large measures of both regular and chaotic orbits is often 
partitioned by such partial obstructions as cantori ||l|] or 
Arnold webs which, although not serving as absolute 
barriers, can significantly impede the motion of a chaotic 
orbit through a connected phase space region. Indeed, 
the fact that, in two-dimensional Hamiltonian systems, 
chaotic orbits can be "stuck" near regular islands for very 
long times was discovered empirically ||3| long before the 
existence of cantori was proven [0 . 

It has also been long known that low amplitude 
stochastic perturbations can accelerate Hamiltonian 
phase space transport by enabling orbits to traverse these 
partial barriers. This was, e.g., explored by Lieberman 
and Lichtenberg [|j, who investigated how motion de- 
scribed by the simplified Ulam version of the Fermi ac- 
celeration map M is impacted by random perturbations, 
allowing for the modified equations (t) 
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(1.1) 



where the "noise" A^ corresponds to a random phase 
shift uniformly sampling an interval [— <y9, +'f]- 



That stochastic perturbations can have such effects on 
Hamiltonian systems is important in understanding the 
limitations of simple models of real systems. In the ab- 
sence of all "perturbations" and any other irregularities, 
the chaotic phase space associated with some idealised 
two- or three-dimensional Hamiltonian system may be 
partitioned into regions which are effectively distinct over 
relatively short time scales. However, even very weak 
perturbations of the idealised model, so small as to seem 
irrelevant on dimensional grounds, can blur these barri- 
ers and permit a single orbit to move from one region to 
another on surprisingly short time scales. 

One practical setting where this may be important is in 
understanding how, in the context of the core- halo model 
g of mismatched charged particle beams, the focusing of 
an accelerator beam can be corrupted by imperfections 
in the magnetic fields. To the extent that such irreg- 
ularities can be modeled as noise, there is the concern 
that noise-induced diffusion can result in particles in the 
beam becoming sufficiently defocused as to hit the walls 
of the container, which is a disaster. Work in this area is 
currently focused on obtaining realistic estimates of the 
noise amplitude and the form of the power spectrum [BJ . 

Another setting is in galactic astronomy. Recent obser- 
vations indicating (i) that many/most galaxies are gen- 
uinely triaxial, i.e., neither spherical nor axisymmetric, 
and (ii) that they contain a pronounced central mass 
concentration suggest strongly that the self-consistently 
determined bulk gravitational potential associated with 



a galaxy contains significant measures of both regular 
and chaotic orbits [|lO|. It was originally expected that, 
in such complex potentials, regular orbits would provide 
the skeleton to support the triaxial structure, and that 
chaotic orbits would serve to fill in the remaining flesh of 
the self-consistent equilibrium |ll[| . However, it appears 
that, in many cases, much of the expected role of regular 
orbits must be played by "sticky" chaotic orbit segments 
since, as a result of resonance overlap, the measure of reg- 
ular orbits in certain critical regions is very small, albeit 
nonzero |l3|. The obvious question is: can low ampli- 
tude perturbations reflecting internal substructures like 
gas clouds and individual stars or the effects of an ex- 
ternal environment destabilise a near-equilibrium on a 
time scale short compared with the age of the Universe? 
Preliminary work would suggest that they can |l3] . 

In both these settings, one knows that weak perturba- 
tions will eventually trigger significant changes in energy 
on some fiducial relaxation time i^, which implies that 
they could have a significant effect. This, however, is not 
the critical issue here. Rather, the question is whether 
low amplitude perturbations can have significant effects 
already on a time scale short compared with the time 
scale on which the value of the energy, or any other iso- 
lating integral, changes significantly. 

In understanding the potential effects of such low am- 
plitude "noise," there are at least three important ques- 
tions which need to be addressed: 

1. How does the effect depend on the amplitude of the 
noise? Is there a threshhold amplitude below which the 
noise is essentially irrelevant, or do the effects turn on 
more gradually? Does the efficacy of the perturbation 
scale as a simple power of the amplitude or does one see 
something more subtle? 

2. To what extent do the details of the noise matter? 
For some problems, such as energy barrier penetration, 
additive (i.e., state- independent) and multiplicative (i.e., 
state-dependent) noises can yield very different results 
||14| . However, the physics here is not the same since 
one is not dealing with a barrier which, in the absence of 
perturbations, is absolute. Rather, one is dealing with an 
entropy harrier ||l^. It would seem that the problem of 
diffusion through cantori or along an Arnold web is more 
similar to problems involving chaotic scattering pq | or 
escapes of unbound orbits from a complex Hamiltonian 



system |17 where, in the absence of perturbations, the 



requisite escape channels exist and it is only a matter of 
how fast any given orbit can find one. 
3. Why does noise lead to accelerated phase space trans- 
port? Granted that the physics is different from diffusion 
in energy, what is the correct physics? One possibility 
is that introducing noise simply fuzzes out the details of 
a purely Hamiltonian evolution that are ensured by Li- 
ouville's Theorem, thus enabling orbits to breach gaps 
which would otherwise be impenetrable. However, some- 
thing very different might be responsible for what is seen. 
This paper aims to address these questions for two- 
dimensional Hamiltonian systems by performing first 



passage time experiments. What this entails is identi- 
fying chaotic orbits which, in the absence of any per- 
turbations, remain "stuck" near regular islands for very 
long times, and determining how the introduction of weak 
noise reduces the escape time. The experiments that were 
performed and interpreted involved both additive and 
multiplicative noise. They also allowed for both white 
noise, which is delta-correlated in time and has a ffat 
power spectrum, and colored noise, which has a finite 
autocorrelation time, so that the power spectrum effec- 
tively cuts off for large frequencies. Finally, the experi- 
ments allowed for both external noise, presumed to exist 
in and of itself, and internal noise, which is accompanied 
by a friction that is related to the noise by a Fluctuation- 
Dissipation Theorem ||l^. To gain additional insights, 
the results of these noisy experiments were also compared 
with experiments in which the unperturbed initial con- 
ditions were evolved in the presence of low amplitude 
periodic driving, so that the breaching of cantori could 
be triggered by modulational diffusion Q. 

Section II describes the experiments that were per- 
formed, and the following three sections report the re- 
sults. Section HI summarises the effects of low amplitude 
periodic driving, indicating the relative importance of the 
amplitude and frequency of the perturbation. Section IV 
describes the effects of different sorts of white noises; and 
Section V generalises the preceding section to the case of 
colored noise. Section VI concludes by summarising the 
evidence that, like periodic driving, noise- induced extrin- 
sic diffusion through cantori is a resonance phenomemon 
which requires substantial power at frequencies compara- 
ble to the natural frequencies of the unperturbed orbit, 
and which has an efficacy that scales logarthmically in 
the amplitude of the perturbation. 



II. A DESCRIPTION OF THE COMPUTATIONS 

The experiments described here were performed for or- 
bits evolved in two representative two-dimensional poten- 
tials, namely the so-called dihedral potential ||l^ for one 
particular set of parameter values, for which the Hamil- 
tonian takes the form 

H = \[pl+pi)-[x'+y') + \{x^+V^f-\x^V^ 

(2.1) 

and the sixth order truncation of the Toda lattice poten- 
tial (20), for which 



h = ^{pI+pI) + 2(^' + y') +^''y-^i 



-x^ + x^y'^ + -y'^ + x'^y + -a;^w^ y^- 

2 ^ 2^ '^ 3 ^ 3^ 



\x'+xV + \x'v'+^^^v'- 



(2.2) 



Extensive explorations of orbits in these Hamiltonians 
would suggest that, in many respects, these potentials 
are generic in the set of nonintegrable potentials admit- 
ting global stochasticity. This is consistent with the fact 
that the experiments performed for this paper yielded 
similar results both potentials. However, these potentials 
are special in the sense that they admit discrete symme- 
tries: the dihedral potential is invariant under a rotation 
by 7r/4; the truncated Toda potential is invariant under a 
rotation by 27r/3. It should be noted for future reference 
that, for relatively low energies, £■ < 40 or so, a charac- 
teristic orbital time scale in each potential corresponds 
to a time i ~ 1 — 3, so that most of the power in typical 
orbits is in frequencies w ~ 1 — 5. 

In both potentials it is easy to find "sticky" [gj chaotic 
orbits which, visually, are very nearly indistinguishable 
from regular orbits for comparatively long times (al- 
though they have short time Lyapunov exponents suffi- 
ciently large that they must be chaotic). Three examples 
are exhibited in the left hand panels of FIG. 1, namely 
two orbits in the dihedral potential, with energies £' = 10 
and E = 20, and an orbit in the truncated Toda potential 
with E = 20. The orbit in FIG.l (a) resembles closely 
what a galactic astronomer would term a regular loop 
orbit; the orbits in FIG. 1 (d) resembles a regular fish. 
The orbit in FIG. 1 (g) is less familiar, but would again 
seem nearly regular. The important point, then, is that if 
the orbit is integrated for a somewhat longer interval, its 
behaviour exhibits an abrupt qualitative change. This 
is illustrated in the center panels of FIG. 1, which ex- 
hibit the same initial conditions, each integrated for an 
interval twice as long. The first two orbits are no longer 
centrophobic, and the third has so changed as to mani- 
fest explicitly the discrete 27r/3 rotation symmetry of the 
truncated Toda potential. In each case, the orbit is far 
more chaotic, as is readily confirmed by the computation 
of a Lyapunov exponent. 

The transition from nearly regular, "sticky" behaviour 
to something more manifestly chaotic occurs once the 
orbit has diffused through one or more cantori that sur- 
round a regular phase space island Q| 0| . The orbit be- 
gins chaotic and remains chaotic throughout, but its ba- 
sic properties exhibit significant qualitative changes after 
the orbit has escaped through the cantori to become "un- 
confined." The precise objective of the work described 
here is to determine how the time required for chaotic or- 
bits to change from sticky to unconfined is altered when 
the orbit is perturbed by low amplitude perturbations. 

Determining the precise location of the outermost can- 
torus is possible, albeit exceedingly tedious [|lj. Fortu- 
nately, however, this is not essential to estimate with 
reasonable accuracy when a "sticky" orbit has become 
"unstuck." Once the orbit has breached the outermost 
confining cantorus, it will typically move quickly to probe 
large portions of the accessible configuration space re- 
gions which were inaccessible before this escape. More- 
over, escape is accompanied by an abrupt increase in 
the value of the largest short time Lyapunov exponent 



p2[ , this reflecting the fact that "sticky" chaotic orbit 
segments confined near regular islands tend to be less 
unstable than unconfined chaotic segments far from any 
regular island |23| . 

As a practical matter, the first escape time for a sticky 
chaotic orbit was identified by (1) using simple polyno- 
mial formulae to delineate approximately the configura- 
tion space region to which the orbit is originally confined, 
and then (2) determining the first time that, with or with- 
out perturbations, the orbit leaves this special region. To 
check that the escape criterion was reasonable, two tests 
were performed: It was verified that, with or without per- 
turbations, small changes in the precise definition have 
only minimal effects on the computed first escape time; 
and that, for the case of unperturbed orbits, the time of 
escape corresponds to a time when the largest short time 
Lyapunov exponent exhibits an abrupt increase. 

This prescription allowed one to identify with reason- 
able accuracy transitions from sticky chaotic to uncon- 
fined chaotic behaviour, but not from chaotic to regular. 
The constant energy surface contains KAM tori, which 
serve as absolute boundaries between regular and chaotic 
behaviour, so that an unperturbed orbit that starts as 
chaotic can never become regular. If, however, the orbit 
is perturbed, the energy is no longer conserved, and it 
becomes possible in some cases for an initially chaotic 
orbit to become regular. 

The experiments described in this paper involved gen- 
erating ensembles of perturbed orbits and then extracting 
statistical properties from these ensembles. In this set- 
ting, two different diagnostics proved especially useful: 

1. The time T(O.Ol) required for one percent of the orbits 
in the ensemble to escape. As described in the follow- 
ing Sections, escapes do not begin immediately. Rather, 
there is typically a relatively extended initial period, the 
duration of which depended on the form of the pertur- 
bation, during which no escapes are observed. Perhaps 
the most obvious number to record would be time when 
the first orbit in the ensemble escaped. However, it was 
found that, in a nonnegligible fraction of the experiments 
that were performed - perhaps 5-10% - one orbit often 
escapes long before any of the others. For this reason, it 
seemed more reasonable to track a diagnostic that is less 
sensitive to comparatively rare exceptions. 

2. The initial escape rate A. In many, albeit not all, 
cases it was found that, once the escape process "turns 
on" at (say) time to i orbits escape in a fashion which, at 
least initially, is consistent with a Poisson process, with 
N{t), the fraction of the orbits which have not escaped, 
decreasing exponentially: 



iV(t) = 7Voexp[-A(t-to)]. 



(2.3) 



The experiments with periodic driving involved solving 
an evolution equation of the form 

—J = - Vy (r) + A sm{ujt + ip)r. (2.4) 



The driving was thus characterised by three parameters, 
namely the frequency w, the ampUtude A, and the phase 
(p. Usually but not always the phase (p was set equal to 
zero. Ensembles of periodically driven orbits were gen- 
erated by (1) specifying a frequency interval [cjjU; + A], 
(2) sampling this interval uniformly to select a collec- 
tion of (usually) 1000 driving frequencies, and then (3) 
integrating the same initial condition with the same am- 
plitude A for each of these frequencies. When looking at 
relatively low frequencies, < w < 100, the frequency 
range was taken to be A = 1.0. For higher frequencies, 
100 < w < 1000, the range A = 10.0. 

The experiments involving intrinsic noise entailed solv- 
ing Langevin equations of the form 






-Vl/(r) - r/v + F, 



(2.5) 



with 77 = ry(v) and F homogeneous Gaussian noise char- 
acterised by its first two moments: 

{Fa{t))^0 and 

{Fa{ti)Fb{t2)}^SabK{v,h-t2), {a,b = x,y). (2.6) 

K{v, t) is the autocorrelation function. For the case of 
delta-correlated white noise. 



if(v,r)=2e77(v),5D(T), 



(2.7) 



where Q denotes a characteristic temperature, the fric- 
tion and noise being related by a Fluctuation-Dissipation 
Theorem. Experiments involving extrinsic noise pro- 
ceeded identically, except that the friction was turned 
off, so that 



'dt^ 



372 = -V^(r) 



(2.8) 



White noise simulations were performed using an al- 
gorithm developed by Griner et al |gj| (see also psf). 
Colored noise simulations were performed using a more 
complex algorithm described in Section V. The experi- 
ments with delta-correlated white noise allowed both for 
additive noise, where 77 is a constant, and multiplicative 
noise, where 77 is a nontrivial function of v. The experi- 
ments with colored noise involved two different choices for 
the form of K{t), in each case allowing for a parameter a 
which characterised the temporal width of the autocor- 
relation function. In every case, ensembles of orbits with 
the same initial conditions were generated by freezing the 
form and amplitude of the noise and performing multiple 
realisations of the same random process using different 
pseudo-random seeds. 



III. PERIODIC DRIVING AND 
MODULATIONAL DIFFUSION 

Experiments involving multiple integrations of the 
same initial condition reveal that escape is (at least) a 



two-stage process: In general there is an initial inter- 
val, often quite extended, during which no escapes occur. 
Only after this interval is there an abrupt onset of escapes 
which, at least for relatively early times, can be well mod- 
eled as a Poisson process, where N{t), the fraction of the 
orbits that have not yet escaped, decreases exponentially. 
An example of this behaviour is illustrated in FIG. 2, 
which was generated for the initial condition exhibited in 
FIG 1 (a), allowing for a frequency interval 2.0 < a; < 3.0 
and an amplitude A = 10~^'^. The straight line exhibits 
a hnear fit to the interval T(O.Ol) < t < 300. 

As asserted already, at time only somewhat larger than 
T(O.Ol), N{t) appears to decrease exponentially. How- 
ever, for t > 400 or so it is clear from FIG. 2 that N{t) 
decreases more slowly. One plausible interpretation of 
this later subexponential decay is that some of the ini- 
tially sticky chaotic orbits have becomes trapped even 
closer to the regular island or, in some cases, have actu- 
ally become regular, so that escape becomes much more 
difficult if not impossible. This interpretation was tested 
by turning off the periodic driving at a late time t = 1024 
and computing both the orbit and an approximation to 
the largest short time Lyapunov exponent for the inter- 
val 1024 < t < 3072. An analysis of the resulting output 
indicated that, in some cases, the orbits which had not 
escaped hy t = 1024 had in fact become regular. 

For fixed frequency interval and phase, both r(O.Ol), a 
measure of the time before escapes begin, and A, the ini- 
tial escape rate once escapes have begun, typically scale 
logarithmically in A, the amplitude of the driving. Six 
examples of this behaviour are provided in FIG. 3, these 
corresponding to the three initial conditions exhibited 
in FIG. 1 for two different frequency intervals, namely 
0.0 < CO < 1.0 and 2.0 < w < 3.0. In each case, the size 
of the error bar has been set equal the difference between 
r(O.Ol) and the time at which the first orbit escapes. In 
most cases, this difference is small, but in some cases it 
becomes appreciable. The fact that the curve is not ex- 
actly linear, and that it levels out for certain ranges of 
amplitude, is not an obvious finite number effect. Dou- 
bling the number of frequencies that were sampled, and 
hence the number of orbits, does not significantly impact 
the overall smoothness of the curve. 

Before the onset of escapes, the rms deviation Sr^ms 
between perturbed and unperturbed orbits typically 



grows as 



Srrms oc ^exp(xt). 



(3.1) 



where A is the driving amplitude and x is comparable to 
the positive short time Lyapunov exponent for the un- 
perturbed orbit. The rms deviation SErms also varies 
linearly with A, but exhibits a much weaker time depen- 
dence. That both these quantities scale linearly in A is 
hardly surprising since periodic driving is a coherent pro- 
cess. The different time dependences reflect the fact that, 
although nearby chaotic orbits tend to diverge exponen- 
tially in configuration space, with or without small per- 
turbations, energy is conserved absolutely in the absence 



of perturbations. That T(O.Ol) scales as log A means that 
escapes begin when Srrms, rather than SErms, assumes 
a roughly constant value, independent of the amplitude 
A. The characteristic value when escapes begin is typ- 
ically Srrms '^ 1 — 2, which implies that the perturbed 
orbits have dispersed to probe most of the region inside 
the confining cantori. 

This leads to a natural interpretation of the escape 
process: Early on, the perturbed orbits remain relatively 
close to the unperturbed orbit, so that it is unlikely that 
they will be able to escape. (The initial conditions were 
so chosen that, in the absence of perturbations, escape 
only occurs at a comparatively late time!) Eventually, 
however, the perturbed orbits will have spread out to 
sample more or less uniformly some region inside the 
bounding cantori. Once this has happened, orbits will 
begin to escape "at random" in a fashion that samples a 
Poisson process. If the holes were very large, one might 
expect that the escape rate at this stage would be nearly 
independent of amplitude. Given, however, that orbits 
still have to "hunt" for tiny escape channels, one might 
expect that A also depends logarithmically on the ampli- 
tude of the perturbation. 

This interpretation is consistent with the expectation 
that an initially localised ensemble of chaotic orbits will 
exhibit an exponential in time approach towards a near- 
invariant distribution that corresponds to a near-uniform 
population of those accessible phase space regions not 
obstructed by cantori pw . It is also qualitatively similar 
to what appears to happen when considering the escape 
of energetically unbound orbits from a complicated two- 
dimensional potential p7| . 

Periodic driving tends to yield the smallest T(O.Ol) 
and largest A for driving frequencies uj comparable to 
the natural frequencies of the unperturbed orbits. For 
example a plot of T(O.Ol) as a function of uj for fixed 
amplitude A and phase ip typically exhibits the smallest 
values of r(O.Ol) for w ~ 1 — 3 and an abrupt increase 
for somewhat larger frequencies. However, low amplitude 
driving can still have an appreciable effect on the time 
of escape even when the driving frequency is much larger 
than the natural frequencies of the unperturbed orbit. 
For example, T(O.Ol) can be significantly shorter than 
the escape time for an unperturbed orbit even for driving 
frequencies as large as w ^ 1000. 

Three examples of how T(O.Ol) varies with cu for fixed 
A and are exhibited in FIG. 4. The three left pan- 
els plot r(O.Ol) as a function oi u; for < uj < 40. The 
three right panels plot r(O.Ol) as a function of log oj for 
I < UJ < 1000. In some cases r(O.Ol) varies smoothly as 
a function of lj for lj ::^ 10; in other cases, considerably 
more irregularity is evident. In either case, however, it is 
apparent that, overall, the efficacy of the driving is set by 
the logarithm of the driving frequency. T(O.Ol) tends to 
increase linearly in log uj. Given the plausible hypothesis 
that this accelerated escape is a resonance phenomenon 
involving a coupling between the driving frequency and 
the natural frequencies of the unperturbed orbits, the 



fact that high frequencies still have an appreciable effect 
can be interpreted as implying that, even though the un- 
perturbed orbit has little power at high frequencies, pe- 
riodic driving can couple via higher order harmonics. 

When, for fixed A and (p, T(O.Ol) and A are compar- 
atively smooth functions of driving frequency, T(O.Ol) 
tends to exhibit only a relatively weak dependence on 
the phase. Different values of tp tend to yield compara- 
ble escape times. If, alternatively, T(O.Ol) depends sen- 
sitively on w, it is more likely that the escape time also 
depends sensitively on ip. However, this trend is not uni- 
form. In some cases, varying ip continuously from to 27r 
changes r(O.Ol) by no more than 10%. In other cases, 
T(O.Ol) can vary by a factor of four, or more. Finally, it 
should be noted that the importance of noise in acceler- 
ating diffusion through cantori can depend sensitively on 
the details of the orbit. Consider, e.g., two initial condi- 
tions in the same potential with the same energy which 
probe nearby phase space regions and which, in the ab- 
sence of perturbations, lead to orbits that escape at com- 
parable times. There is no guarantee that ensembles of 
periodically driven orbits generated from these different 
initial conditions and evolved with the same amplitudes, 
phases, and driving frequencies will exhibit similar val- 
ues of r(O.Ol) and A, even if the unperturbed orbits have 
power spectra that are almost identical. 



IV. WHITE NOISE 

For stationary Gaussian white noise with zero mean, 
everything is characterised by the autocorrelation func- 
tion K{ti — t2) — 077(r, v)(5£)(ii — ^2), the form of which 
is determined in turn by 77. Choosing 77 to be constant 
yields additive noise. Allowing for a nontrivial depen- 
dence on r or V yields multiplicative noise. One aim of 
the work described here was to determine the extent to 
which the detailed form of the noise matters. This was 
done by first performing experiments involving additive 
noise, and then comparing the results with experiments 
that involved multiplicative noise of two forms, namely 
77 (X w^ and 77 ex v^"^ , where v denotes the orbital speed. 
The importance of friction was tested by comparing ex- 
periments that included a friction related to the noise by 
a Fluctuation-Dissipation Theorem ||l^ with experiments 
with no friction at all. 

If the friction and noise are to mimic internal degrees 
of freedom that are ignored in a mean field description, 
one anticipates on dimensional grounds that the temper- 
ature O will be comparable to a typical orbital energy. 
For this reason, most of the experiments that were per- 
formed, including those described here, involved freez- 
ing the temperature at a value Q ^ E and exploring the 
effects of varying the amplitude of 77. The relative nor- 
malisations of the multiplicative and additive noises were 
fixed by setting 



^(v) =iio{v/{v)) 



±2 



(4.1) 



where 770 denotes the constant rj appropriate for addi- 
tive noise and (v) denotes the average speed of the un- 
perturbed orbit. Comparing additive and multiphcative 
noise entailed comparing experiments with the same 770. 

Considering only two forms of multiplicative noise in- 
volves probing the tip of an iceberg: other multiplica- 
tive noises could in principle have very different effects. 
However, the two cases examined here do allow one to 
ask whether the overall effect of the friction and noise 
can change significantly if one allows the statistics of the 
noise to vary along an orbit. The particular forms cho- 
sen here were motivated by two considerations: (1) If the 
noise is intended to mimic discreteness effects in a plasma 
or a galaxy (i.e., electrostatic or gravitational Rutherford 
scattering), the friction should depend on velocity P8| . 
(2) Allowing for a relatively strong dependence on speed, 
ry oc w^^, should make even relatively small differences 
comparatively easy to see. 

Overall, the effects of white noise are very similar to 
the effects of periodic driving. In particular, escape was 
again observed to be a two stage process, involving an 
initial interval during which different realisations of the 
same initial condition diverge inside the confining can- 
tori, followed by an abrupt onset of escapes which, at 
least initially, is well approximated as a Poisson process. 
Moreover, as for the case of periodic driving, N{t) de- 
creases subexponentially at late times, possibly because 
some of the noisy orbits have become regular or, at least, 
more closely trapped near a regular island. 

FIGURE 5 exhibits plots of log N{t) generated for one 
representative initial condition, corresponding to the or- 
bit in FIG. 1 (a). Here the experiments involved additive 
noise and friction related by a Fluctuation-Dissipation 
Theorem, with a fixed 6 = 10 and IQ-^ < ryo < 10"''. It 
is evident that, as rjQ decreases and the friction and noise 
become weaker, the escape time T(O.Ol) increases and 
the escape rate A decreases. 

More careful examination reveals further that, for fixed 
0, both r(O.Ol) and A scale logarithmically in tjq. This 
is illustrated in FIG. 6, which exhibits T(O.Ol) and a best 
fit value of A for the ensembles used to construct FIG. 5 
(along with some ensembles with intermediate values of 
77). This logarithmic dependence can be understood by 
analogy with what was observed for periodic driving if 
one notes that, in the presence of noise, the rms devia- 
tion between perturbed and unperturbed orbits typically 
scales as 



Srrms (X (776)^/2 exp(xt), 



(4.2) 



a conclusion that can be motivated theoretically |29| and 
and has been confirmed computationally Q. 

It is interesting that this two-stage evolution - an 
epoch without escapes followed by an epoch with es- 
capes apparently sampling a Poisson process ~ can also 
be observed in the absence of noise if one considers a 
strongly localised ensemble of initial conditions trapped 
near a regular island and ascertains the time at which 



each member of the ensemble escapes. For example, 
an ensemble of orbits sampling a cell of size 0.002 cen- 
tered about the initial condition used to generate FIG. 6 
yielded T(O.Ol) = 310 and A = 0.000587, which should 
be compared with the values r(O.Ol) = 131 and A = 
0.00153 resulting for a single initial condition evolved 
with ry = IQ-^. 

Perhaps the most significant conclusion about white 
noise is that, at least for the examples considered here, 
the details are largely irrelevant. For fixed 8 and 770, 
the values of the escape time T(O.Ol) and the decay rate 
A arc both essentially the same for the simulations with 
additive noise and those with multiplicative noise with 
77 oc v^^. The computed values of T(O.Ol) and A are also 
nearly independent of whether or not one allows for a 
friction term. The only significant differences between 
simulations with and without friction arise at late times 
when the energies of individual orbits have changed ap- 
preciably from their initial values. In this case, allowing 
for a friction term to counterbalance the noise assures 
that, overall, the energies of the orbits exhibit smaller 
changes, so that the ensembles evolved with noise tend 
to have somewhat smaller changes in energy. 

An example of this insensitivity is provided in FIG. 7, 
which was generated from orbits with the initial condi- 
tion of FIG. 1 (a), with 6 = 10 and 770 = 10"^ The solid, 
dot-dashed, and triple-dot-dashed curves exhibit N{t) for 
ensembles evolved in the presence of both friction and 
noise, incorporating, respectively, additive noise, multi- 
plicative noise with 770 oc u^ , and multiplicative noise with 
770 oc u~^. The dashed curve corresponds to an ensemble 
evolved with additive noise in the absence of friction. The 
obvious point is that, for a very long time, these curves 
are nearly indistinguishable. 

This insensitivity is again consistent with the hypoth- 
esis that noise-induced diffusion through cantori is a res- 
onance phenomenon, and that the only thing that mat- 
ters is that the noise have significant power at frequencies 
comparable to the natural frequencies of the unperturbed 
orbit. If, however, one alters the form of the autocorre- 
lation function so as to suppress power at frequencies 
comparable to these natural frequencies, one would ex- 
pect that the effects of the noise should decrease, so that 
r(O.Ol) increases and A decreases. The extent to which 
this is true is discussed in Section V. 



V. COLORED NOISE 

The objective of the experiments described here was to 
explore the effects of random perturbations with autocor- 
relation times sufficiently long that they cannot be mod- 
eled as delta-correlated white noise. Once again it was 
assumed that the noise is stationary and Gaussian with 
zero mean and, for simplicity, attention was restricted to 
noise that is additive. However, it was no longer assumed 
that K{t) is delta-correlated in time. 



Attention focused on two types of colored noise: The 
first is generated by the Ornstein-Uhlenbeck process (see, 
e.g., jl^l), and is characterised by an autocorrelation 
function K{t) that decays exponentially, i.e.. 



^W = "^loi exp(-a|T|). 



(5.1) 



The second involves an exponential decay modulated by 
power law: 



K 



(r) = a^oiexp{^a\T\){l + a|r| + ^r^) . (5.2) 



The corresponding spectral densities are, respectively. 



S{u;) = -^ 



and 



S{co) 



8^2 



Stt (cj2+a2)3' 



(5.3) 



(5.4) 



The autocorrelation times are tc — l/a and tc — 2/a. 
In both cases, white noise corresponds to a singular limit 
with a — >• oo and a'^^i -^ cx3, but il^i/a -^ const. As 
for the case of multiplicative noise, these two examples 
only probe the tip of an iceberg. However, an analysis 
of their effects does provides insight into the question of 
how a finite autocorrelation time can impact phase space 
transport in a complex phase space. 

Generating white noise numerically is comparatively 
straightforward, requiring little more than producing a 
sequence of pseudo-random impulses. Generating colored 
noise takes more thought. The algorithm exploited here 
was motivated by the recognition that, in the context of 
a stochastic differential equation, a white noise random 
process X{t) can serve as a source to define a colored 
noise process Y{t). As a concrete example, consider how 
Gaussian white noise can be used to implement a random 
process with an autocorrelation function given by (5.2). 

Given one stochastic process, X(t), one can define a 
second stochastic process, Y{t), implicitly as a solution 
to the stochastic differential equation 



d^Y{t) 



3a 



d'^Yjt) 
dt^ 



3a- 



■ dYjt) 
dt 



a^Y{t) 



X{t). 
(5.5) 



Since the coefficients in eq. (5.5) are time-independent 
constants and X{t) is stationary, it is clear that, if this 
equation be solved as an initial value problem, at suffi- 
ciently late times Y{t) can also be considered stationary 
provided only, as is true, that the dynamical system (5.5) 
is stable. 

By expressing X(t) and Y(t) in terms of their Fourier 
transforms, it is easy to see that, neglecting the effects of 
nontrivial boundary conditions (e.g., choosing boundary 
conditions at to -^ — oo), the spectral densities Sx{^) 
and Syi^) for the two processes satisfy 



Sy{^) 



Sx{uj) 



(u^ 



a 



2^3 ■ 



(5.6) 



Assuming, however, that X(t) corresponds to white 
noise, Sx{i^) = <^x i^ ^ constant, so that the stochas- 
tic process with spectral density Sy necessarily corre- 
sponds to colored noise. Indeed, by performing an in- 
verse Fourier transform it becomes evident that 5y(w) 
corresponds to the stochastic process (5.2) with 



37rcr|/(8a5). 



(5.7) 



A colored random process defined by the Langevin equa- 
tion (2.5) or (2.8), with an autocorrelation function of the 
form (5.2), is equivalent mathematically to a collection 
of white noise processes. Solving (5.5) for Y(t) yields the 
colored input required to solve (2.5) or (2.8). 

The one remaining question involves normalisations. 
To compare different colored noises with each other or 
with an appropriately defined white noise limit, one must 
decide what should be meant by noise with variable au- 
tocorrelation time but fixed amplitude. This was done 
here by considering sequences of random processes with 
different values of a and, for each a, selecting a^^i such 
that 



K^{T)dT 



K,. 



'hitt 



■iT)dT. 



(5.8) 



In other words, fixed amplitude but variable autocorrela- 
tion time was assumed to correspond to different colored 
noises for which the time integral of the autocorrelation 
function assumes the same value. Noting that 



(5.9) 



Kwhite{T)dT = 297?, 

it follows that, for the stochastic process (5.2), 
3a77e ' -2 



Kir) 



■exp(-a|r|)(l-Ka|r| + yr2). (5.10) 



The Ornstein-Uhlenbeck process requires a normalisation 
K{t) = m]Q exp(-a|T|). (5.11) 

In the experiments modeling intrinsic noise, the col- 
ored noise was augmented by a friction ry, and the fric- 
tion and noise were related by a linear Fluctuation- 
Dissipation Theorem in terms of a temperature Q. In the 
experiments modeling extrinsic noise, the friction van- 
ished and O had no independent meaning as a tempera- 
ture. Aside from a, which fixes the autocorrelation time, 
all that matters is the quantity rjQ, which sets the am- 
plitude of the noise. 

As for white noise, the evolution of an ensemble of 
noisy colored orbits is a two stage process. After an ini- 
tial epoch without escapes, during which different mem- 
bers of the ensemble diverge exponentially, escapes turn 
on abruptly, with the first percent of the orbits escaping 



within an interval T(O.Ol) much shorter than the time 
before the first escape. Interestingly, though, the sec- 
ond phase is often more complex than what is observed 
for white noise. Instead of evolving in a fashion that is 
well fit pointwise by an exponential decrease, the number 
remaining, N{t) often exhibits "plateaux" and "jumps" 
with (almost) no, and especially large, decreases. (A de- 
tailed examination of the data reveals that such irregu- 
larities can also arise for white noise and periodic driving, 
but that are usually much less conspicuous in that case.) 

FIGS. 8-11 exhibit N{t) generated for an initial con- 
dition where such irregularities are comparatively small. 
FIGS. 8 and 9 exhibit data generated for different 4800 
orbit ensembles generated, respectively, for the stochastic 
processes (5.1) and (5.2), for the same initial condition 
as FIG. 5, evolved with 8 = 10.0, a = 0.2, and variable 
7]. FIGS. 10 and 11 exhibit analogous plots for the same 
initial condition with 77 = 10~^ and variable a. FIG. 12 
exhibits data for a second initial condition in the same 
potential with the same energy for which the early-times 
irregularities are especially conspicuous. 

These irregularities imply that a pure exponential fit 
is often not justified. Moreover, even when such a fit 
is justified, one finds that, for fixed a, significantly dif- 
ferent values of 77 can yield very similar slopes [ pO| . In 
this sense, it is not accurate to state unambiguously that 
the escape rate scales logarithmically with amplitude. 
However, what does remain true is that, overall, escapes 
tend to happen more slowly in the presence of lower am- 
plitude perturbations, and that any systematic ampli- 
tude dependence is very weak, certainly much weaker 
than a simple power law oc ri~P with p of order unity. 
Moreover, even though the observed escape rates exhibit 
considerable irregularities, the one percent escape time 
r(O.Ol) does not. As for the case of white noise and 
periodic driving, r(O.Ol) scales logarithmically in ampli- 
tude. Several examples are exhibited in FIGS. 13 (a) 
and (c), which plot T(O.Ol) as a function of log 77 for two 
different initial conditions. In each case, the diamonds 
represent an Ornstein-Uhlenbeck process and the trian- 
gles the stochastic process (5.2). In panel (a) a — 2.0; in 
panel (c) a = 0.2. 

As for white noise, one also finds that, seemingly in- 
dependent of a, the presence or absence of friction is 
largely irrelevant. This is, e.g., evident from Table 1 
which, for two values of a, namely a = 00 (white noise) 
and a — 0.02 (autocorrelation time r = 50), exhibits 
T(O.Ol) as a function of log?] both in the presence and 
the absence of friction. 

Another obvious conclusion is that the efficacy of col- 
ored noise is a decreasing function of a, the quantity that 
sets the autocorrelation time t^. When a is very large, 
so that tc is extremely short, color has virtually no ef- 
fect. However, as a decreases and ic increases, T(O.Ol) 
increases. In particular, for a <C 1, this corresponding to 
an autocorrelation time that is long compared to a char- 
acteristic orbital time scale, r(O.Ol) is typically much 
longer than what is found in the white noise limit. This 



behaviour is evident from FIGS. 13 (b) and (d) which 
exhibit T(O.Ol) as a function of log a for fixed r] — Vir'^ . 
As for FIGS. 13 (a) and (c), the diamonds and triangles 
represent, respectively, the stochastic processes (5.1) and 
(5.2). The horizontal dashed line represents the white 
noise value towards which the data converge for a ^ 00. 
The obvious inference from this, and other, plots is that, 
the dependence of r(O.Ol) on a or t^ is again roughly 
logarithmic. This is reminiscent of the fact that, as dis- 
cussed in Section HI, the efficacy of periodic driving tends 
to scale logarithmically in the driving frequency. 

Determining the overall efficacy of colored noise as 
a source of accelerated phase space transport thus in- 
volves an interplay between amplitude and autocorrela- 
tion time, each of which, in the "interesting" regions of 
parameter space contributes logarithmically to r(O.Ol) 
and (modulo the aforementioned caveats) A. 



VI. DISCUSSION 

Just as for diffusion triggered by low amplitude peri- 
odic driving, the overall efficacy of noise-induced diffu- 
sion of "sticky" chaotic orbits scales logarithmically in 
the amplitude of the perturbation. For both white and 
colored noise, the one percent escape time T(O.Ol) scales 
logarithmically in the amplitude of the perturbation and, 
at least for white noise, so does the initial escape rate. 
The details of the perturbation seem largely unimpor- 
tant: The presence or absence of a friction term appears 
immaterial, and allowing for (at least some forms of) mul- 
tiplicative noise also has comparatively minimal effects. 
For the case of white noise the only thing that seems to 
matter is the amplitude. 

This logarithmic dependence on amplitude implies that 
what regulates the overall efficacy of the noise in inducing 
phase space transport is how fast perturbed noisy orbits 
diverge from the original unperturbed orbit. Escapes en- 
tail a two stage process, namely (i) an early interval dur- 
ing which noisy orbits diverge from the unperturbed orbit 
without breaching cantori and (ii) a later interval during 
which, in many cases, orbits escape seemingly at random 
in a fashion that samples a Poisson process. That T(O.Ol) 
scales logarithmically in amplitude reflects the fact that 
escapes typically begin once Srrms, the rms separation 
between perturbed and unperturbed orbits, approaches a 
critical value comparable to the size of the region in which 
the "sticky" orbits are originally stuck. That the escape 
rate tends in many cases to exhibit at least a rough loga- 
rithmic dependence reflects the fact that, even after the 
noisy orbits have spread out to sample a near-invariant 
population inside the confining cantori, random kicks can 
facilitate phase space transport by helping the orbits to 
"find" holes in the cantori. 

That escapes begin when 5rrms, rather than SErms, 
approaches a critical value has an important implica- 
tion for how one ought to envision the escape process. 



Naively, it is not completely obvious whether friction and 
noise facilitate phase space transport by "jiggling" indi- 
vidual orbits or, since they change the orbital energy, by 
"jiggling" the effective phase space hypersurface in which 
the orbits move. The fact that Srmis sets the scale on 
which things happen demonstrates that, in point of fact, 
the former interpretation is more natural. 

Allowing for colored noise can significantly reduce the 
rate of phase space transport, but only when the au- 
tocorrelation time tc becomes comparable to, or larger 
than, a characteristic crossing time. If the noise is such 
that the spectral density S{ijj) has only minimal power 
at frequencies comparable to the natural frequencies of 
the unperturbed orbit, its efficacy in inducing accelerated 
phase space transport will be reduced significantly. For 
a sufficiently large and tc sufficiently small, the two col- 
ored noises that were explored yield essentially the same 
results as did white noise; and similarly, for a sufficiently 
small and t^ sufficiently long, the effects of the noise must 
become essentially negligible, so that one recovers the be- 
haviour observed for an unperturbed orbit. For interme- 
diate values, however, the escape statistics do depend on 
the value of a. Moreover, this dependence is reminiscent 
of the effects of periodic driving in at least one impor- 
tant respect: For the case of periodic driving, quantities 
like T(O.Ol) exhibit a roughly logarithmic dependence on 
the driving frequency to. For the case of colored noise, 
T(O.Ol) exhibits a roughly logarithmic dependence on a. 

This suggests strongly that, like modulational diffusion 
triggered by periodic driving, noise-induced phase space 
diffusion is intrinsically a resonance phenomenon. If the 
Fourier transform of the noise has appreciable power in 
the frequency range where the unperturbed orbit has ap- 
preciable power, noise-induced diffusion will be compar- 
atively efficient. If, however, the noise has little power 
at such frequencies, it will serve as a much less efficient 
agent for phase space transport, although the effects need 
not be completely negligible. 

To the extent that, as suggested by the numerical 
experiments described here, the details are relatively 
unimportant, the effects of noise as a source of phase 
space transport are determined by two physical quanti- 
ties, namely (i) the amplitude and (ii) the autocorrelation 
time. Increasing the amplitude makes noise more impor- 
tant; increasing the autocorrelation time makes noise less 
important. 

As a concrete example, consider stars orbiting in an 
elliptical galaxy comparable in size to the Milky Way 
but located in the central part of a cluster like Coma, 
where the typical distance between galaxies is only five 
to ten times larger than the diameter of a typical galaxy. 
Here there are two obvious sources of noise which one 
might consider, namely (i) "discreteness effects" reflect- 
ing the fact that the galaxy is made of individual stars 
rather than a dustlike continuuum and (ii) the near- 
random influences of the surrounding environment. Dis- 
creteness effects result in gravitational Rutherford scat- 
tering, which can be modeled reasonably [E§] by friction 



and delta-correlated white noise related by a Fluctuation- 
Dissipation Theorem where, in natural units, 77 ~ 10^^ — 
10^^. To the extent that the near random influences of 
the surrounding environment can be attributed primarily 
to a small number of neighbouring galaxies, it is also easy 
to estimate their amplitude and typical autocorrelation 
time. Given that the nearest neighbouring galaxy is typi- 
cally separated by a distance ~ 5 — 10 times the diameter 
of the galaxy in question, and that the relative velocities 
of different galaxies in a cluster are usually comparable to 
the typical velocities of stars within an individual galaxy, 
one expects that the autocorrelation time i* is of order 
5—10 characteristic orbital times tcr- Presuming, how- 
ever, that the perturbing influences of nearby galaxies 
reflect tidal effects, their overall strength should scale as 
D~^ , where D is the distance from the galaxy in question, 
so that a typical amplitude 77 ^ 10^'^ — 10^^. 

In this setting, discreteness effects give rise to com- 
paratively weak noise with a very short autocorrelation 
time. Environmental effects given rise to a considerably 
stronger noise with a much longer autocorrelation time. 
The longer autocorrelation time tends to suppress the ef- 
fects of environmental noise, but, even so, it would seem 
likely that, as a source of accelerated phase space trans- 
port, environmental noise will be significantly more im- 
portant than discreteness noise. 
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FIG. 1. (a) A chaotic initial condition with E = 10 evolved in the dihedral potential for a time t = 512. (b) The same 
orbit integrated for t = 1024. (c) The power spectra |a::(<^)| and |j/(i^)| for the orbit in (a), (d) A chaotic initial condition with 
E — 20 evolved in the dihedral potential for a time t — 512. (e) The same orbit integrated for t — 1024. (f) |a;(i.<j)| and 12/(1^)1 
for the orbit in (d). (g) A chaotic initial condition with E = 20 evolved in the truncated Toda potential for a time t — 300. 
(h) The same orbit integrated for t = 600. (i) |x(a;)| and \y{u>)\ for the orbit in (g). 
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FIG. 2. N{t), the fraction of the orbits from a 4001 orbit 
ensemble not yet having escaped at time t, computed for the 
initial condition exhibited in Fig. 1 (a), allowing for a per- 
turbation of amphtude A — 10^^"' with variable frequencies 
2.0 < Lu < 3.0. 
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FIG. 3. (a) T(O.Ol), the first escape time for 1% of an 
ensemble of 1000 integrations of the initial condition of Fig. 
1 (a) , driven with frequencies < tJ < 1 , plotted as a function 
of the logarithm of the amplitude A of the perturbation, (b) 
The same for an ensemble with 3 < c^; < 4. (c) T(O.Ol), the 
first escape time for 1% of an ensemble of 1000 integrations 
of the initial condition of Fig. 1 (d) , driven with frequencies 
< w < 1, plotted as a function of the logarithm of the am- 
plitude A of the perturbation, (d) The same for an ensemble 
with 3 < a; < 4. (e) r(O.Ol), the first escape time for 1% of an 
ensemble of 1000 integrations of the initial condition of Fig. 1 
(g), driven with frequencies < a; < 1, plotted as a function 
of the logarithm of the amplitude A of the perturbation, (f) 
The same for an ensemble with 3 < cj < 4. 
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FIG. 4. (a) r(O.Ol), the first escape time for 1% of an 
ensemble of 1000 integrations of the initial condition of Fig. 
1 (a), driven with amplitude A — 10^^'^, plotted as a function 
of frequency range u; < Q < u> + 1 for < cj < 40. (b) The 
same information for 1 < w < 1000, now plotted as a function 
of logtj. The dashed line represents the escape time for the 
unperturbed orbit, (c) T(O.Ol), the first escape time for 1% 
of an ensemble of 1000 integrations of the initial condition of 
Fig. 1 (c), driven with amplitude A — 10^^'^, plotted as a 
function of frequency range ui < Q, < u> + 1 for < w < 40. 
(d) The same information for 1 < uj < 1000, now plotted as 
a function of logcj. The dashed line represents the escape 
time for the unperturbed orbit, (e) r(O.Ol), the first escape 
time for 1% of an ensemble of 1000 integrations of the initial 
condition of Fig. 1 (e), driven with amplitude A — 10^^'^, 
plotted as a function of frequency range ui < fl < cj + 1 for 
< u> < 40. (f) The same information for 1 < oj < 1000, now 
plotted as a function of log ui. The dashed line represents the 
escape time for the unperturbed orbit. 
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FIG. 5. N{t), the fraction of the orbits from a 2000 or- 
bit ensemble not yet having escaped at time t, computed 
for the initial condition exhibited in Fig. 1 (a), allowing for 
additive white noise with Q = 10 and variable rj = lO^'* 
(broad dashes), rj = 10~^ (triple-dot-dashed), r; — 10^® 



(dot-dashed), r; = 10 ^ (narrow dashes), rj 
and r? = lO"'* (dots). 
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FIG. 6. (a) r(O.Ol), the first escape time for 1% of an en- 
semble of 2000 white noise integrations of the initial condition 
of Fig. 1 (a), with = 10 and variable rj. (b) A, the rate 
at which orbits in this ensemble escape, fitted to the interval 
r(O.Ol) < t < 256. 
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FIG. 7. N{t), the fraction of the orbits from a 2000 orbit 
ensemble not yet having escaped at time t, computed for the 
initial condition exhibited in Fig. 1 (a) with O = 10 and 
770 = 10~^. The four curves represent additive white noise 
and constant rj (solid), additive white noise with no friction 
(dashed), multiplicative noise with rj (x v^ (dot-dashed), and 
multiplicative noise with rj oc v^'^ (triple-dot-dashed). 



FIG. 9. N{t), the fraction of the orbits from a 4800 orbit 
ensemble not yet having escaped at time t, computed for the 
initial condition exhibited in Fig. 1 (a), allowing for friction 
and colored noise given by eq. (5.2) with a — 0.2, O — 10, and 
variable rj = 10"" (dashed), ri = 10"^ (dot-dashed), 77 = 10"*^ 
(double-dot-dashed), 77 = 10~^ (brad-dashed), and rj = 10^* 
(solid). 
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FIG. 8. N{t), the fraction of the orbits from a 4800 orbit 
ensemble not yet having escaped at time t, computed for the 
initial condition exhibited in Fig. 1 (a), allowing for friction 
and additive Ornstein-Uhlenbeck noise with a — 0.2, O = 10, 
and variable r; = 10^" (dashed), rj = 10~^ (dot-dashed), 
rj = 10^® (double-dot-dashed), rj — 10~^ (broad-dashed), and 
r]= 10"* (solid). 



FIG. 10. N{t), the fraction of the orbits from a 4800 orbit 
ensemble not yet having escaped at time t, computed for the 
initial condition exhibited in Fig. 1 (a), allowing for friction 
and Ornstein-Uhlenbeck noise, with O = 10, 77 = 10"'', and 
either white noise (dot-dashed) or variable q = 20 (dashed), 
a = 2.0 (triple-dot-dashed), a = 0.2 (broad dashes), and 
a = 0.02 (solid) 
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FIG. 11. N{t), the fraction of the orbits from a 4800 or- 
bit ensemble not yet having escaped at time t, computed for 
the initial condition exhibited in Fig. 1 (a), allowing for 
friction and colored noise given by eq. (5.2) with O — 10, 
7] = 10~^, and either white noise (dashes) or variable a = 20 
(dot-dashed), a = 2.0 (triple-dot-dashed), a — 0.2 (broad 
dashes), and a = 0.02 (solid) 
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FIG. 13. (a) r(O.Ol), the first escape time for 1% of an 
ensemble of 4800 integrations, computed for the initial condi- 
tion used to generate FIG. 1 (a), plotted as a function of log?) 
for fixed a = 2.0 for the stochastic processes defined by (5.1) 
(diamonds) and (5.2) (triangles), allowing for both friction 
and noise, (b) r(O.Ol) for the same initial condition, plotted 
as a function of logo for fixed rj — 10^"", for the stochastic 
processes (5.1) (diamonds) and (5.2) (triangles), again allow- 
ing for both friction and noise. The dashed line represents the 
asymptotic value for white noise (a —> oo). (c) The same as 
(a), albeit for a different initial condition and with a — 0.2. 
(d) The same as (b) , albeit for the initial condition in (c) and 
with T] = 10^ . as (a) and (b) for another initial condition. 



TABLE I. The 1% escape time T(O.Ol) for orbits in an en- 
semble evolved in the dihedral potential (2.1) with the initial 
condition used to generate Fig. 1 (a), setting B = 10.0 and 
allowing for variable rj. The colored noise was generated by 
the stochastic process (5.2). 



FIG. 12. The same as FIG. 8 for a different initial condi- 
tion. 



log r] a = oo a = oo a = 0.02 a = 0.02 

with friction no friction with friction no friction 



4 


21.621 


21.672 


71.930 


71.811 


5 


41.740 


41.660 


113.024 


112.980 


6 


71.863 


71.816 


172.396 


172.393 


7 


82.744 


82.724 


277.684 


277.682 


8 


131.368 


131.616 


277.704 


277.709 



